use dataset_cces, clear

eststo clear
eststo: quietly reg  envindex  anywells i.year, vce(cluster zipcode)
estadd local ZFE = "No"
eststo: quietly areg envindex anywells i.year, a(zipcode) vce(cluster zipcode)
estadd local ZFE = "Yes"
eststo: quietly areg envindex anywells ideoq i.year, a(zipcode) ///
	vce(cluster zipcode)
estadd local ZFE = "Yes"
eststo: quietly areg envindex wellsideo2 wellsideo3 wellsideo4 ///
	wellsideo5 anywells q2 q3 q4 q5 i.year, a(zipcode) cl(zipcode)
estadd local ZFE = "Yes"

label variable anywells "Any wells"
label variable ideoq "Liberalism"
label variable wellsideo2 "Any wells X liberalism = 2"
label variable wellsideo3 "Any wells X liberalism = 3"
label variable wellsideo4 "Any wells X liberalism = 4"
label variable wellsideo5 "Any wells X liberalism = 5"
label variable q2 "Liberalism = 2"
label variable q3 "Liberalism = 3"
label variable q4 "Liberalism = 4"
label variable q5 "Liberalism = 5"

#delimit;

esttab 
	using "_output/tableA15.tex"
	,
		label
		replace
		booktabs
		nomtitles
		collabels(none)
		drop(*year*) 
		cells(b(star fmt(3)) se(par))
		stats(ZFE N, label("Zip code fixed effects" "Observations") 
			fmt(%20.0fc))
		nonote
		;

#delimit cr

